Astronomy & Astrophysics manuscript no. s09 


©ESO 2010 


July 30, 2010 





Role of emission angular directionality in spin determination of 
accreting black holes with a broad iron line 



J. Svoboda, M. Dovciak, R. Goosmann* and V. Karas 



o 



OS 
(N 

< 

o 

Oh- 

6 

CO 

> 

OO 

m 

OO 

o 

On 
O 



X 



Astronomical Institute, Academy of Sciences, Bocm II 1401, CZ-14131 Prague, Czech Republic 
Received 24 February 2009 / Accepted 4 August 2009 

ABSTRACT 

Aims. The spin of an accreting black hole can be determined by spectroscopy of the emission and absorption features produced in 
the inner regions of an accretion disc. In this work, we discuss the method employing the relativistic line profiles of iron in the X-ray 
domain, where the emergent spectrum is blurred by general relativistic effects. 

Methods. Precision of the spectra fitting procedure could be compromised by inappropriate accounting for the angular distribution of 
the disc emission. Often a unique profile is assumed, invariable over the entire range of radii in the disc and energy in the spectral band. 
An isotropic distribution or a particular limb-darkening law have been frequently set, although some radiation transfer computations 
exhibit an emission excess towards the grazing angles (i.e., the limb brightening). By assuming a rotating black hole in the centre 
of an accretion disc, we perform radiation transfer computations of an X-ray irradiated disc atmosphere (NOAR code) to determine 
the directionality of outgoing X-rays in the 2-10 keV energy band. Based on these computations, we produce a new extension to 
the KY software package for X-ray spectra fitting of relativistic accretion discs. 

Results. We study how sensitive the spin determination is to the assumptions about the intrinsic angular distribution of the emitted 
photons. The uncertainty of the directional emission distribution translates to ^ 20% uncertainty in the determination of the marginally 
stable orbit. We implemented the simulation results as a new extension to the KY software package for X-ray spectra fitting of 
relativistic accretion disc models. Although the parameter space is rather complex, leading to a rich variety of possible outcomes, 
we find that on average the isotropic directionality reproduces our model data to the best precision. Our results also suggest that 
an improper use of limb darkening can partly mimic a steeper profile of radial emissivity. We demonstrate these results in the case 
of XMM-Newton observation of the Seyfert galaxy MCG-6-30-15, for which we construct confidence levels of^ 2 statistics, and on 
the simulated data for the future X-ray IXO mission. Our simulations, with the tentative IXO response, show a significant improvement 
that can qualitatively enhance the accuracy of spin determination. 

Key words. Accretion, accretion-discs - Black hole physics - Galaxies: active - X-rays: binaries 



1. Introduction 

There is now strong evidence suggesting that the radiation of ac- 
tive galactic nuclei as well as some compact binary stars in the 
Galaxy originates in part from th e gas of an accretion disc or - 
biting around a central black hole dReeslll998 l :lKato et alX l998). 
A great deal of in formation about these objects has been ob- 
tained via X-rays (iPetersonll 1 997t lKrolikl[T999h . Therefore, reli- 
able models are needed to describe production of high-energy 
photons emerging from the accretion disc surface and subse- 
quent propagation of these photons towards a distant observer. 

Directional distribution of the outgoing radiation is among 
the important aspects that must be addressed. Limb darkening 
traditionally refers to the gradual diminution of intensity in the 
image of the surface of a star as one moves from the centre of 
the image to the edge. It is a consequence of uneven angular dis- 
tribution of the radiation flux eme rging from the stellar surface 
dChan drasekha rl 1 960t iMihalasI 1 9781) . Limb darkening results as 
a combination of two effects: (i) the density of the surface lay- 
ers decreases in the outward direction, and (ii) temperature also 
drops as the distance from the centre of the star increases. The 
outgoing radiance is therefore distributed in a non-uniform man- 
ner. The actual form of the angular distribution depends on the 
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physical mechanism responsible for the emission and the geo- 
metrical proportions of the source. The limb darkening law is 
widely applied also to describe radiation coming from an ac- 
cretion disc around a black hole. In this case, relativistic effects 
can significantly enhance the observed anisotropy by bending 
the light rays and boosting the photon energy in the direction of 
motion of the emitter. As a result of this interplay between local 
physics of light emission and the global effects of the gravita- 
tional field, different processes contribute to the final (observed) 
directional anisotropy of the emission. 

We consider the angular distribution of the emitted X-rays 
in the context of accreting black holes. The source is an equa- 
torial accretion disc in which the energy is drawn by con- 
version of gravitational energy o f the orbiting matter gradu- 
ally sinking i nto the black hol e dShakura & Svunvaevl ll973: 
Pringle 1981: [Frank etal.1 12002), and electromagnetic process 
creati ng coronal flares and illuminating the underlying accretion 
disc (Gal eev et al.lll979t iNavaksMrj 12000: ICzernv et alJl200H 
and further references cited therein). The relevant objects are 
accreting supermassive black holes in nuclei of galaxies and 
stellar-mass black holes in compact binary systems (see, e.g., 
iMcClintock & Remillardl2006tlMillerl2007l:lNandra et al.l2007l 
for broad reviews of the subject). The form of the X-ray spec- 
trum and, in particular, the discrete features at 6-7 keV provide 
a powerful tool to explore the physical properties of the emitting 
region, close to the black hole horizon. The line profile often 
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appears to be smeared and asymmetric, with the red wing ex- 
tending down to ~ 3 keV or even less. Spectral modelling can 
help us to determine the system parameters, nam ely, the black 
hole spin dRevnolds et dJl2005t Miller et alJl2009h . 

It is important to realise that the emission angular direction- 
ality probes the physical conditions of the emitting medium. 
Apart from the above-mentioned temperature stratification, it 
depends also on the opacity profile. The latter is typically a 
frequency-dependent quantity, hence the observed radiation also 
forms at a variable physical depth within the disc. Therefore, the 
limb darkening/brightening property depends on the wavelength. 
We take this into account by maintaining the energy resolution 
in our computations. 

The resulting radiation is dominated by X-rays in a few kilo- 
electronvolt (keV) band, which come from the innermost regions 
of the accreting system and bear imprints of a st rong gravita- 
tional field and rapi d orbital motion at their origin (Fabia n et alj 
2000; Fabian 2008). The spectrum is a combination of a multi- 
colou r thermal component from the accr eting medium of the 
disc dPage&Thornelll974t iLi et al.ll2005h . the power-law com- 
ponent originating from the disc corona dGaleevetal.lll979t 
ICzernv et al.l 12004). and the reflection and absorption features 
arising on th e disc surface and in the extended corona alon g the 
line of sight dFabian et alJI 19891: iRevnolds & Nowakll2003l) . For 
active galaxies, the thermal multi-colour black body component 
dominates in the UV and rarely reaches the sub keV energies. 
This component, however, occurs at higher energy of the spec- 
trum of black hole binaries. 

This paper is organised as follows. Section [2] formulates the 
model and presents its basic assumptions. Section [3] explores 
how the constraints on the black hole spin depend on the direc- 
tionality of the spectrum emerging around the iron line energy. 
To this end we construct the confidence level diagrams of the x 2 
statistics which demonstrate the level of expected sensitivity of 
the model. First we explore simple analytical approximations of 
the limb-darkening profile, which do not depend on energy, and 
then we consider numerical results of our self-consistent radia- 
tive transfer computations of both the continuum and the spec- 
tral features between ~ 2-10 keV. We demonstrate the results by 
"fitting" the artificially simulated data, for which all parameters 
are perfectly under our control, and by reanalysing the XMM- 
Newton observation of the Seyfert 1 galaxy MCG-6-30-15. We 
discuss our results in section [4] and briefly summarise them in 
section [5] 

2. Model assumptions and requisites 

2.1. Gravitational field near a rotating black hole 

According to general relativity, the gravitational field of a ro- 
tating black hole influences the intrinsic spectrum of a nearby 
source by ch anging the energy of photons and bending the shape 
of light rays (Misner et al] |1973l) . The assumed geometry of the 
source, i.e. the accretion disc with high orbital velocity near its 
inner rim, is another factor contributing to significant anisotropy 
of the emission in the observer's frame. We remark that geomet- 
rical optics provides an appropriate framework describing these 
effects near astrophysical black holes, where the wavelength of 
X-rays is much shorter than the typical length-scale of the sys- 
tem represented by the black hole gravitational radiusQ In the in- 
nermost regions of accretion discs, bulk motion of the gas takes 

1 The gravitational radius is r g = GM/c 2 « 1.48 x 10 13 Mg cm, where 
the central black hole mass is expressed in terms of M 8 =M/(10 8 Mq). 
The velocity of the Keplerian orbital motion near a black hole is v K ~ 



place at a significant fraction of the speed of light, indicating 
that special relativistic aberration and Doppler boosting should 
be important. 

The gravitation al field is described by the Kerr metric 
dMisneretalJll973L chapt. 33): 

, AS , A sin 2 6 i \2 S , , 

ds 2 = - — df 2 + — - — (dcp-Lodt) + -dr 2 + Ed6» 2 , (1) 

where the Boyer-Lindquist spheroidal coordinates (f, r, 9, (f>) and 
geometrised units (c = G = 1) are used. Metric functions 
are: A(r) = r 2 - 2r + a 2 , £(r, 6) = r 2 + a 2 cos 2 9, A{r, 9) = 
(r 2 + a 2 ) 2 - A(r) a 2 sin 2 9, and oj(r, 0) = 2ar/A(r, 9); a denotes the 
specific rotational angular momentum (spin) of the central body. 
These functions are assumed not to be perturbed by gravitational 
effects of the external matter (accretion disc self-gravity is ne- 
glected). Such an assumption is well established if we consider 
the inner regions of the disc, which within the distance of several 
gravitational radii contai ns only very low mass compared to the 
central black hole mass dGood man 2003; Ka ras et aUl2004l) . 
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Fig. 1. Orbital velocity *V(r) of co-rotational motion near a ro- 
tating black hole, as given by formula ([3]) for three values of 
the black-hole dimensionless angular momentum parameter a 
(spin). The thick part of each curve indicates the range of radii 
above the marginally stable orbit, r > r ms (a), where the circular 
motion is stable. The thin curve indicates an unstable region at 
small radii. 



Unlike weakly gravitating, non-compact bodies (e.g. main- 
sequence stars), black holes strongly affect light passing near 
their horizon. This influence bears not only the signature of black 
hole mass but it also depends on its rotation via frame-dragging. 
Although the latter effect is rather weak, it does shape the details 
of the final spectrum and hence can be employed to measure 
the black hole angular momentum. The conversion factor from 
the angular momentum 7 p h ys (in physical units) to the angular 
momentum J (in geometrical units) reads: J = (G/c 3 )7 p hy S . The 
geometrised dimension of J is the square of the length [cm 2 ]. 

It is convenient to make all geometrised quantities dimen- 
sionless by scaling them with the appropriate power of mass 

2.1 x 10'°(r/r g ) _1/2 cm s~', and the corresponding orbital period is T K » 
3.1 x 10 3 (r/r g ) 3/2 Mg s (here we neglect an order-of-unity correction 
factor due to the relativistic frame-dragging effect which, however, will 
be included later in the paper). 
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M. The dimensionless specific angular momentum, a = J/M 2 , 
spans the range — 1 < a < 1, where the positive/negative value 
refers to the motion co/counter-rotating with respect to the black 
hole. We will further assume co-rotational motion only (a > 0). 
We can scale all lengths by M to reach dimensionless units. 

For the Keplerian angul ar velocity of the orbital motion, we 
obtain (Bardee n et al.l[l97l 



1 



.3/2 



+ a 



(2) 



For the linear velocity with respect to a locally non-rotating ob- 
server, we have (see Figure [T]) 



<V(f) = 



r 2 - 2ar 1/2 + a 2 

ai/2 Jym + a ) ■ 



(3) 



The velocity at the marginally stable orbit reaches a consid- 
erable fraction of the light speed c and has a similar value 
'V(r) m 0.5 - 0.6 c for any value of the angular momentum. For 
large spin, a small dip develops in the velocity p rofile near the 
horiz on. Although it is an interesting feature (see Stuchlf k et al.l 
l2005h . its magnitude is far too small to be recognised in our anal- 
ysis of the complicated structure of the model x 2 space. 

Notice that the above-mentioned interval of a-values follows 
from the assumption, traditionally accepted by astronomers, that 
we are dealing with a Kerr black hole, for which the outer hori- 
zon is located at r+ = 1 + (1 — a 2 ) 1 ^ 2 . Hence, the magnitude of a 
is thought to be less than unity in order to have a regular horizon 
and avoid the case of naked singularity, although the latter pos- 
sibility cannot be straightforwardly rejected just on the basis of 
spectra fitting. 

The standard disc scenario assumes that the accretion disc 
stops at the innermost stable circular orbit (ISC O, also called 
the m arginally stable orbit), r m = r ms (a), where (B ardeen et alJ 

[1971 



= 3+Z 2 -[(3-Z 1 )(3+Z 1 +2Z 2 )]\ 



(4) 



Zi = l+(l-fl 2 )5[(l+a)5+(l-a)5] andZ 2 = (3a 2 +Z 2 )^. Notice 
that r ms (a) spans the range of radii from r ms = 1 for a = 1 (the 
case of maximally co-rotating black hole) to r ms — 6 for a = 
(static black hole). As discussed further below, Laor's (Il99ll) 
model adopts the standard disc scheme and further assumes that 
rotation of the black hole is limited by an equilibriu m value, 
q=0.9 98, because of capture of photons from the disc dThornd 
119741) . This implies r ms =1.23. 



2.2. Black-hole signatures in accretion disc radiation 

The physical parameters of the black hole - namely, its an- 
gular momentum a - are imprinted in the spacetime geom- 
etry and they influence the velocity of the orbiting mate- 
rial as well as the radiation propagating through that space- 
time. Thus, the observed spectrum is affected - the equiv- 
alent width and the centroid energy of spectral lines differ 
from the v alue expected on the basis of a purely Newtonian 
model (e.g. Martocchia et al. 2000; Dabrowski & Lasenbv 2001; 
iMiniutti et al.ll2003HVaughan & Fabianll2004 and further refer- 
ences cited therein). In other words, the relativistic effects of- 
fer the possibility to reveal the presence of a compact body in 
the centre of the accretion disc, and to determine its parameters. 
This assertion is based on an implicit assumption that we can 
constrain the local physics forming the intrinsic spectrum of the 
source with sufficient accuracy. 



The effects of strong gravity may be concealed from us if 
the intrinsic spectrum of the accretion disc is not known well 
enough. This is particularly an issue with respect to the angu- 
lar momentum determination because the geometrical effect of 
frame-dragging due to a rotating black hole is a rather subtle 
one. In order to assess potential inaccuracies, we investigate how 
sensitive the determination of the a parameter is to the angular 
distribution of the disc emission, i.e., the limb-darkening (or the 
limb-brightening) law in the spectral range of 2 < E < 10 keV. 

Numerical as well as semi-analytical methods were devel- 
oped and combined to study the spectral features from relativis- 
tic accretion discs, especially the profiles of the iron-line com- 
plex between 6-7 keV. An important practical attribute from the 
view point of applicability of any particular method i s its im- 
pleme ntation within the XSPEC data analysis package (lArnaudl 
1996). In this context, over almost two decades the most widely 
used mode l of the relativistic disc spectral line is the one by 
lLaon (ll99n. which includes the effects of a rapidly rotating Kerr 
black hole. The LAOR model sets the dimensionless angular mo- 
mentum a to the canonical value of a — 0.998, so it cannot 
be the subject of a r igorous data fitting procedure (although see 
ISvoboda et al]l2008h . 

Although the mathematical properties of the extreme Kerr 
spacetime (a 2 = 1) are significantly different from the sub- 
extreme case (a 2 < 1), the graphs of the redshift factor as well 
as the emission angle are practically identical. Therefore we do 
not expect any detectable differences of the observed spectra be- 
tween these two cases. Secondly, the main effect on the standard 
accretion disc spectrum and the relativistic broad line formation 
arises from the assumed link betwee n the inner edge of the dis c 
and the black hole s pin parameter dRevnolds & Fabian 2008). 
Svoboda et a D (I2008h have examined the accuracy of the spin 
fitting using the LAOR model. They demonstrate that, although 
such an approach is not self-consistent, the accuracy of the re- 
sultin g a-values is actually quite good. 

iDovciak et al.1 d2004bl) have relaxed the limitation of the 
LAOR model regarding the fixed value of the spin. These au- 
thors allowed for a to be fitted to the dat a by y 2 minimisa- 
tion in a suite of KY models. Further more. iBeckwith & Done 
(2004) have developed the KDLINE, and Brenneman & Reynolds 
(2006) have developed the KERRDISK models which are en- 
dowed with similar functionality. The latter authors performed 
useful tests demonstrating that KY and KERRDISK give com- 
patible results when they are set to equivalent parameter val- 
ues. Numerical codes have been developed i n dependently by 
several other groups (Martocchia et a □ 120001; [Viergut3 11993b 
Zakh arov&Repinll2004t iFuerst & Wul 12004 ICadez & Calvanl 



2005h using different techniques. The LAOR kernel has been re- 
cently applied by IMiniutti et a l. (2007) who reanalysed the pro- 
totypical source MCG-6-30-15 using Suzaku observations to 
determine the relativistic blurring and to measure the black hole 
para meters. Their res ults confirm the high value of the spin. 

iMiller et"aT] (l2009h applied the modelling of relativistic spec- 
tral features to estimate the black hole spin from a sample of 
stellar-mass black holes. They used the joint constraints from 
both the disc reflection and continuum and found evidence for a 
broad range of black hole spin parameters. 

The angular emissivity law, Ai(/i e , r e , E e ), defines the distri- 
bution of the intrinsic intensity outgoing from each radius r e of 
the disc surface with respect to the perpendicular direction. The 
emission angle e = arccos fi e is measured from the disc normal 
direction to the equatorial plane, in the disc co-moving frame, 
i.e. in the local Keplerian frame orbiting with the angular ve- 
locity Q K (r). Likewise the intrinsic energy E e is measured with 
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respect to the local frame. The total disc emission can be written 
in the form of product 



I(r e ,p. e ,E e ) = , R(r & )fi e M(iu e ,r e ,E e )8{E e ), 



(5) 



where the radial part is well approximated by a power law, 

■R(r e ) = rl q (q = const), (6) 

and &(E e ) is the energy dependence of the locally emitted ra- 
diation. When needed we will also use the broken power-law 
profile which is a combination of two profiles (|6]l matched at the 
transition radius, also coded in the KY set of models. 
The redshift factor g and the emission angle 6> e are 



C 



arccos 



(7) 



where S= l+ar^ 3 ^ 2 , C = l-3r~ l +2ar~ 3 t 2 ;t; and rj are constants 
of motion connected with symmetries of the Kerr space time. 

The observed radiation flux is then obtained by integrating 
the intrinsic emission over the entire disc surface, from the in- 
ner edge (r = r m ) to the outer edge (r - r out ), weighted by the 
transfer function T(r e , <f> e , O , a) determining the impact of rela- 
tivistic energy change (Doppler and gravitational) as well as the 
lensing effect for a distant observer directed along the inclina- 
tion angle 6^ (see|Cunninghamlll975l : lAsaokalfl989t iKaras et all 
ll992tlKarasll2006l) . 

Given the high velocity of the orbital motion and the strong- 
gravity light bending near the black hole, the effect of directional 
anisotropy of the local emission is enhanced. For this reason it is 
important to describe the angular distribution in a correct man- 
ner; ad hoc choices of the limb-darkening law may lead to errors 
in the determination of the model best fit parameters, including 
the inaccuracy in a parameter which are difficult to control, or 
they may prevent us from estimating the statistical confidence of 
the model. In the case of black hole accretion discs, this compli- 
cation becomes important because the aberration, beaming and 
light-bending effects grow rapidly towards the inner edge of the 
disc. Therefore, even a small discrepancy between the assumed 
and the correct angular emissivity profiles becomes greatly en- 
hanced in the observer frame. 



2.3. Effects of the emission angular directionality 

In optics, Lambert's cosine law describes an emitter producing a 
radiation intensity that is directly proportional to yt/ e . Lambertian 
surfaces exhibit the same apparent radiance when viewed from 
any angle O . Likewise, Lambertian scattering refers to the situ- 
ation when the surface radiates as a result of external irradiation 
by a primary source and the scattered light is distributed accord- 
ing to the same cosine law. This is, however, a very special cir- 
cumstance; directionality of the emergent light is sensitive to the 
details of the radiation mechanism. For example, the classical 
result of the Eddington approximation for stellar atmospheres 
states that the effective optical depth of the continuum is t = |, 
and so the emergent intensity is described by the limb-darkening 

law, 7(jU e ) oc yU e + §. 

In the case of a fluorescence iron line produced by an il- 
luminated plane-parallel slab, t he angular d i stribution was in- 
vestigated by various authors (Bas koi 119781: lG eoree & Fabian 
[T99lIlHaardllT993tlGhisemni et alJl994tlRevnolds et aTB oOO). 
In that case a complicated interplay arises among the angular 
distribution of the primary irradiation, reflection and scattering 
in the disc atmosphere. Several authors pointed out that it is es- 
sential for the reliable determination of the model parameters 



to determine t he angular direct i onality of the broad line emis- 
sion correctly. Martocch ia et al.l (12000) noted, by employing the 
lamp-post model, that "...the broadening of the observed spectral 
features is particularly evident when strongly anisotropic emis- 
sivity laws, resulting from small h [i.e., the lamp-post elevation 
above the equatorial plane], are considered." 

The important rol e of the emission angular d irectionality was 
clearly spelled out bv lBeckwith & Pond d2004l) : "...the angular 
emissivity law (limb darkening or brightening) can make sig- 
nificant changes to the derived line profiles wh ere light bend- 
ing is im porta nt" (see their Fig . 9-13). Similarly, JDpvciak et al. 
( 2004ajb|) and iBeckwifhl (120051) compared the relativistic broad 
lines produced under different assumptions about the emission 
angular directionality. However, to verify the real sensitivity of 
the models to the mentioned effect of directionality, it is nec- 
essary to connect the radiative transfer computations with the 
spectral fitting procedure, and to carry out a systematic analysis 
of the resulting spectra, taking into account both the line and the 
continuum in the full relativistic regime. Here we report on our 
resul ts from such compu tations. 

iRevnolds et al.l (120041 sec. 4.3) argue that the combined ef- 
fect of photoelectric absorption in the disc and Compton scat- 
tering in the corona more affect the iron line photons emerg- 
ing along grazing light rays than continuum photons. They con- 
clude that the line equivalent width should be diminished for 
observers viewing the accretion disc at high inclination an- 
gle s. Such a tr end c an be seen also in the l amp-p ost model 
of iMatt et al l (|1992|) and iMartocchia eTal] (|2000, see their 
Figure 11). However, in the latter work this diminution is less 
pronounced when we compare it with the case of intrinsically 
isotropic emissivity. 

More recently, Niedzwiecki & Zvcki (2008) studied the ef- 
fect of different limb-darkening laws on the iron-line profiles. 
They point out that the role of emission directionality can be 
quite significant once the radial emissivity of the line is fixed 
with sufficient confidence. However, this is a serious assump- 
tion. In reality, the radial emissivity is not well constrained by 
current models. 

The angular dependence of the outgoing radiation is deter- 
mined by the whole interconnection of various effects. We de- 
scribe them in more detail below (sec. 13. 4t . Briefly, the conclu- 
sion is such that realistic models require numerical computations 
of the full radiative transfer. We have developed a complete and 
consistent approach to such radiative transfer computations in 
the context of the broad iron-line modelling together with the un- 
derlying continuum computations. As described in considerable 
detail below, we performed the extensive computations which 
are necessary in order to reliably determine the impact of the 
emission angular anisotropy on spectral fitting results (namely, 
on the determination of the black hole angular momentum). In 
particular, we describe detailed results from the investigation of 
the model goodness (by employing an adequate statistical anal- 
ysis of the complicated^ 2 parameter space). Such the analysis 
has not been performed so far in previous papers because de- 
tailed stepping through the parameter space and proper re-fitting 
of the model parameters was not possible due to the enormous 
complexity of the models and extensive computational costs. 

In regular stars and their accretion discs, the relativistic 
effects hardly affect the emerging radiation. The situation is 
very different in the inner regions of a black hole accretion 
disc, where the energy shift and gravitational lensing are sig- 
nificant. The observed signal can be boosted or diminished 
by the Doppler effect combined with gravitational redshift: 
I(fi e ) / I(fi ) = g 3 , where the g values span more than a decade. 
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Fig. 2. Contours of redshift factor g(r,4>) near a maximally rotating black hole, a = 1, depicted in the equatorial plane (x,y). The 
black hole and the accretion disc rotate counter clock-wise. The inner region is shown up to r — 20 gravitational radii from the black 
hole (denoted by dark circle of unit radius around the centre). A distant observer is located towards the top of the figure. Two cases 
of different observer inclinations are shown. Left: O - 30 deg. Right: <9 = 70 deg. The colour bar encodes the range acquired by 
g(r, <p), where g > 1 corresponds to blueshift (approaching side of the disc), while g < 1 is for redshift. 



15 30 45 60 75 90 15 30 45 60 75 90 




Fig. 3. Contours of the local emission angle, 6 e (r, <p). The same parameters as in the previous figure. The local emission angle spans 
the entire range, from to 90 degrees. This is due to the combined effects of aberration and light bending which grow greatly near 
the inner rim of the disc. In the inner region the photons are boosted in the direction of rotation and they emerge along grazing 
angles. Asymptotically, 9 e (r, (f>) — > O for r — > oo. 
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Fig. 4. Left: directional distribution of the intrinsic emissivity (see sec. l3.1l for details). Right: an example of the directional distri- 
bution from the numerical computations, showing results from the reprocessing model with the continuum photon index r = 1 .9, 
integrated over 2-100 keV energy range (solid line), and in several different energy sub-ranges (line styles are indicated in the 
inset). The latter graph demonstrates the presence of limb-brightening effect in comparison with respect to the isotropic emission. 
This effect is clearly visible as a function of the emission angle fi e , although its magnitude is smaller than the limb-brightening 
approximation used in the left panel. See sec. l3.4l for details. 



Figure|2]shows the redshift factor g near the black hole, proving 
that general relativistic effects indeed change the photon energy 
significantly. 

As men tioned above, many authors hav e adopted the defin- 
ing choice ( Chandrasek harl[T960t lLaorll 19911) of the cosine pro- 
file for the line angular emissivity: Ai(fi e , r e , E e ) — 1 + 2.06 fi e . 
This relation describes the energy-independent limb-darkening 
type of profile. However, the choice is somewhat arbitrary in the 
sense that the physical assumptions behind this law are not satis- 
fied at every radius over the entire surface of the accretion disc. It 
has been argued that the limb-darkening characteristics need to 
be modified, or even replaced by some kind of limb brightening 
in the case of X-ray irradiated disc atmospheres with Compton 
reflec tion dHaardtll 19931 iGhisellini etaDll994t IZvcki & Czernvl 
11994 . The latter should include the energy dependence, as the 
Compton reprocessing of the reflected component plays a signif- 
icant role. The angle-dependent computations of the Compton 
reflec tion demonstrate these effects convincingly (ICzernv et al.l 
120041) . Indeed, the same effect is seen also in our computations, 
as shown in the right panel of Figure |4] The increase of emis- 
sivity with the emission angle strongly depends on the ionisa- 
tion state of the reflecting material, so the actual situation can be 
quite complex (Goosmann et al. 2007b. 

A question arises of whether the current determinations of 
the black hole angular momentum might be affected by the un- 
certainty in the actual emissivity angular distribution, and to 
what degree. In fact, this may be of critical concern when future 
high-resolution data become available from the new generation 
of detectors. We have therefore carried out a systematic investi- 
gation using the KY code to reveal how sensitive the constraints 
on the dimensionless a parameter are with respect to the possible 
variations in the angular part Ai(/J. & ) of the emissivity. 

Figure [3] shows the contours of the local emission angle 6 e 
from the accretion disc, taking into account their distortion by 
the central black hole and assuming that the emitted photons 
reach a distant observer at a given view angle 8 . Although the 
contours are most dramatically distorted near the horizon, where 
both aberration and light bending effects grow significantly, the 



emission angle is visibly different from the observer inclination 
even quite far from the horizon, at a distance of several tens r g . 
This is mainly due to special-relativistic aberration which decays 
slowly with the distance as the disc obeys Keplerian rotation at 
all radii. 

In summary, figures [2]-[3] demonstrate the main attributes of 
the photon propagation in black hole spacetime relevant to our 
problem: the energy shift and the direction of emission depend 
on the view angle of the observer as well as on the angular mo- 
mentum of the black hole. Notice that, near the inner rim, the 
local emission angle is indeed highly inclined towards the equa- 
torial plane where at the same time the outgoing radiation is 
boosted. These effects are further enhanced by gravitational fo- 
cusing, which we also take into account in our calculations. 

It should be noted that the energy shifts and emission an- 
gles near a black hole, such as those shown in Figs. 2-3, have 
been studied by a number of authors in mutually complemen- 
tary ways. The figures shown here have been produced by plot- 
ting directly the content of FITS format files t hat are encode d 
in the corresponding KY XSPEC routines (see iDovcia kl (12004 . 
where an atlas of contour plots is presented for different in- 
clinations and spins). A nalogous figures were also shown in 
iBeckwith & Dona (12004 . who depict the dependence of the co- 
sine of the emission angle on the energy shift of the received 
photon and the emission radius. 

The above quoted papers concentrated mainly on the discus- 
sion of the energy shifts and the emission directions of the in- 
dividual photons, or they isolated the role of relativistic effects 
on the predicted shape of the spectral line profile. They clearly 
demonstrated that the impact of relativistic effects can be very 
significant. However, what is still lacking is a more systematic 
analysis which would reveal how these effects, when integrated 
over the entire source, influence the results of spectra fitting. To 
this end one needs to perform an extensive analysis of the model 
spectra including the continuum component, match the predicted 
spectra to the data by appropriate spectra fitting procedures, and 
to investigate the robustness of the fit by varying the model pa- 
rameters and exploring the confidence contours. 
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One might anticipate the directional effects of the local emis- 
sion to be quite unimportant. The argument for such an expecta- 
tion suggests that the role of directionality should grow with the 
source inclination, whereas the unobscured Seyfert 1 type AGNs 
(where the relativistically broadened and skewed iron line is usu- 
ally expected) are thought to have only small or moderate incli- 
nations. However, this qualitative trend cannot be used to quanti- 
tatively constrain the model parameters and perform any kind of 
precise analysis, needed to determine the black hole spins from 
current and future high-quality data. Such an analysis has not 
been performed so far, and we embark on it here for the first 
time. 

3. Iron Ka line band examined with different 
directionalities of the intrinsic emissivity 

3.1. Approximations to the angular emission profile 

We describe the methodology which we adopted in order to ex- 
plore the effects of the spectral line emission directionality. To 
this end we first employ simple approximations, neglecting any 
dependence on the photon energy and the emission radius. We 
set the line intrinsic emissivity from the planar disc to be de- 
scribed by one of the following angular profiles, 



Case 1: 



Case 2: M(p e ) = 
Case 3: 



ln(l+^ 1 ) (Haardtl993) 
1 (locally isotropic emission) 
1 + 2.06 yu e (Laor 19911 



(8) 



and 



with the radial profiles being a power law 7?(r e ) = 
£(£e) = <5(£ e - Eq). The three cases correspond, respectively, to 
the limb-brightened, isotropic, and limb-darkened angular pro- 
files of the line emission. 

The limb-brightening law by HaardU d 19931) describes the an- 
gular distribution of a fluorescent iron line emerging from an 
accretion disc that is irradiated by an extended X-ray source. 
The relation was obtained from geometrical considerations 
and agrees well with more detailed Mon te-Carlo computations 
dGeorge & Fabianl [l991: Matt et al.ll 19921) . The physical circum- 
stances relevant for the limb-darkening law are different, and we 
include this case mainly because it is implemented in the laor 
model and frequently used in the data analysis. The isotropic 
case, dividing all limb-brightening and limb-darkening emissiv- 
ity laws, is included in our analysis for comparison. Notice that 
an extra geometrical factor // e is involved due to simple geomet- 
rical projection of the disc surface in eq. (0. 

The radial profile of the emission is set to a unique power 
law, eq. ©, over the entire range of radii across the disc. The 
directionality formula ([8]) of the intrinsic emissivity and the re- 
sulting spectral profiles are illustrated in Fig. [4] (left panel). 
Naturally, more elaborate and accurate approximations have 
been discussed in the literature for some time. For example, 
Ghise llini etalJ(fl99l in their eq. (2) include higher-order terms 
in fi e to describe the X-ray reprocessing in the single-scattering 
Rayleigh approximation. However, at this stage the first-order 
terms are sufficient for us to demonstrate the differences between 
the three cases. Later on we will proceed towards numerical radi- 
ation transfer computations that are necessary to derive realistic 
profiles of the emission angular distribution and to keep their 
energy dependence. 

The dependence of the spectral profiles on the angular distri- 
bution A1(// e ) of the intrinsic emission is shown in figures [5]-[6] It 
is apparent from Figure [5] that the limb brightening case makes 



the line profile broader (in the left panel) and the height of the 
blue peak lower (in the right panel) than the limb darkening case 
for the same set of parameters, which can consequently lead to 
discrepant evaluation of the spin. While Figure [5] shows the line 
profile for an extended disc, Figure [6] deals with a narrow ring. 
In this case, a typical double-horn profile develops. Although 
the energy of the peaks is almost entirely insensitive to the emis- 
sion angular directionality, the peak heights are influenced by the 
adopted limb darkening law. In the right panel of Figure [6] the 
ratio of the two peaks is shown for the three cases of the emis- 
sion angular directionality according to eq. ([8]). The influence of 
the directionality is apparent and it is comparable to the effect 
of the inclination angle. The spin value has only little impact for 
large inclinations. 

In general, we notice that the extension and prominence 
of the red wing of the line are indeed related to the intrin- 
sic emission directional ity. This was examined in detail by 
Beck wifh & Pond (120041) . who plotted the expected broad-line 
profiles for different angular emissivity and explored how the 
red wing of the line changes as a result of this undetermined 
angular distribution. Similarly, the inferred spin of the black 
hole must depend on the assumed profile to a certain degree. 
This is especially so if the spin is determined by the extremal 
redshift of the line's red w ing. It has been recently argued 
dNiedzwiecki & Zvc ki 2008) that the directionality effect could 
be a real issue for the accuracy of line fitting if other parameters 
(such as the radial emissivity profile) are well constrained. 

However, simple arguments are insufficient to assess how in- 
accurate the spin determination might be in realistic situations, 
as the spectral fitting procedure involves several components ex- 
tending over a range of energy above and below the iron-line 
band. We illustrate this in Figure [7j noticing that the different 
prescriptions produce very similar results outside the line en- 
ergy, but they do differ at the broad line energy range. The the- 
oretical (background-subtracted) profiles of the relativistic line 
cannot alone be used to make any firm conclusion about the er- 
ror of the best fit parameters that could result from the poorly 
known angular emissivity. To this end one has to study a con- 
sistent model of the full spectrum. With a real observation, the 
sensitivity to the problem of directionality (as well as any other 
uncertainty inherent in the theoretical model) will depend also 
on the achieved resolution, energy binning and the error bars of 
the data. 



3.2. Example: MCG-6-30-15 reanalysis 

We reanalysed a long XMM-New ton observation of a nearby 
Seyfert 1 galaxy MCG-6-30-15 dFabian et alj l2002h to test 
whether the different directionality approximations can be dis- 
tinguished in the current data. The observation took place in 
summer 2001 and the acquired exposure time was about 350 ks. 
The spectra reveal the presence of a broad and skewed iron line 
(Fabi an et ail l2002t iBaUantvne et alj 120031: IVaughan & Fabian! 
2004). We reduced the EPN data from three sequential revo- 
lutions (301, 302, 303) using the SAS software version 7.1.20 
We used standard tools for preparing and fitting the data con- 
tained in the HEASOFT software package version 6.40 Using 
the FTOOL MATHPHA (part of the HEASOFT) we joined the 
three spectra to improve the statistical significance. Further, we 
used the XSPEC version 12.2 for the spectral analysis. We re- 
binned all the data channels to oversample the instrumental en- 



2 http://xrnrn.esac.esa.int/sas 

3 http://heasarc.gsfc.nasa.gov 
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Fig. 5. Left: theoretical profiles of the relativistic line (the KYRLINE model, without continuum), corresponding to the three cases 
in the left panel of Fig. [4] The lines are normalised with respect to the height of the blue peak. Model parameters are a = 0.9982, 
q — 3, r m = r ms (a) = 1.23, r out = 400, o = 30 deg, Eq = 6.4 keV. Right: the same as in the left panel, but with the normalisation set 
in such a way that the radiation flux is identical in all three profiles. 
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Fig. 6. Left: The same as in Fig. [5] but for a narrow ring with r m = 4.7, r out = 4.8 (other parameters are a = 0.9982, q — 3, 
o = 30 deg, Eo = 6.4 keV). The lines are normalised with respect to the height of the blue peak. Right: The ratio of the two peak 
heights for the three cases of the emission angular directionality according to eq. © and for the two values of the spin and the 
inclination angle (crosses and circles are for O = 30 deg; empty and filled triangles are for 8 = 60 deg). 



ergy resolution maximally by a factor of 3 and to have at least 
20 counts per bin. The first condition is much stronger with 
respect to the total flux of the source - 4 x 10~ n ergcnT 2 s _1 
(1.1 x 10 6 cts) in the 2-10keV energy interval. 

We appl i ed the same continuum model as presented in 
Fabi an et all (120021) : the power law component (photon index 
r = 1.9) absorbed by neutral hydrogen (column density «h = 
0.41 x 10 21 cm -2 ). This simple model is sufficient to fit the data 
above as 2.5 keV, which is also satisfactory for our goal of re- 
producing the overall shape of the broad iron line. Other com- 
ponents need to be added to the model in order to fully un- 
derstand the spectrum formation and to decide between viable 
alternatives, including the presence of outfl ows and a combi- 
nation of absorption a nd reflection effects dMiller et al.| [2008; 
iTurner &Millerl 12009). However, our goal here is not to give 
precedence to any of the particular schemes. Instead, we rely on 
the model of a broad line and we test and compare different an- 
gular laws for the emission. 



The residuals from the simple power law model can be ex- 
plained by a complex of a broad line and two narrow lines. Due 
to the complexity of the model, it is hard to distinguish between 
the narrow absorption line at E » 6.75 k eV and an emission 
line at E * 6.97 keV dFabian et al.ll2002l) . Although the level 
of x 2 values stays almost the same, the parameters of the broad 
iron line do depend on the continuum model and the presence 
of narrow lines. By adding an absorption line at £ a; 6.75 keV, 
the best fit rest energy of the broad component comes out to 
be E — 6.7 keV, and so the broad line origina tes in a moder- 
ately i onised disc. This result is consistent with Balla ntvne et alJ 
(120031) . 

More details of our reanalys i s of th e broad iron line shape 
are described in I Svoboda et alJ (120081) . Here, we present the 
results of the one-dimensional steppar command in Figure [8] 
which demonstrates the expected confidence of a-parameter 
best-fit values for the two extreme cases of directionality, 
Case 1 and Case 3. The model used in the XSPEC syntax is: 

PHABS*(POWERLAW+ZGAUSS+ZGAUSS+KYRLINE). The fixed 
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Fig. 7. Left: unfolded spectrum generated by the POWERLAW + KYRLINE model with parameters of the line identical as in previous 
figure. We assumed the limb-brightening profile, i.e. Case 1 of eq. ®. Right: the ratio of the model on the left to the same model 
with the Case 3 (limb-darkening) profile. The normalisation of the line was allowed to vary during the fit. The changed from 
1.04 to 1.30 between the two cases. The details about the simulation of the data points are the same as described in the section [331 




0.8 0.85 0.9 0.95 1 0.8 0.85 0.9 0.95 1 

spin a/M spin a/M 

Fig. 8. The best-fit values of;f 2 statistics for the s pin parameter, wh ich we obtained by gradually stepping a from 0.8 to 1. XMM- 
Newton data for MCG-6-30-15 were employed (Fa bian et alJl2002T) . Left: the limb-brightening profile (Case 1). Right: the limb 
darkening profile (Case 3). The dashed line is the 90% confidence level (2 sigma). See the text for a detailed description of the 
model. 



parameters of the model are the column density «h = 0.41 X 
10 21 cm 2 , the photon index of power law T = 1.9, the red- 
shift factor z = 0.008, the energy of the narrow emission line 
E em = 6.4 keV, and the energy of the narrow absorption line 
£abs = 6.77 keV. The parameters of the broad iron line were 
allowed to vary during the fitting procedure. Their default val- 
ues were £broad = 6.7 keV for the energy of the broad iron line, 
9 = 30deg for the emission angle, q\ = 4.5, qi — 2 and 
r b = 10r g for the radial dependence of the emissivity (the ra- 
dial part of the intensity needs to be rather complicated to fit the 
data and can be expressed as a broken power law: fi(r e ) = r e c " 
for r e < r b , and 7?(r e ) = r e 92 for r e > rb). 

The determined best-fit values for the spin are virtually the 
same for both cases, independent of the details of the limb- 
brightening/darkening profile. However, this result arises on ac- 
count of the growing complexity of the model. The differences 
between the two cases become hidden in different values of the 



other parameters - especially in q\, qi and rb, i.e. the parame- 
ters characterising the radial dependence of the line emissivity in 
KYRLINE as a broken power law with a break radius r b . We find: 
(i) Abroad = 6.60(1), o = 31.5(7)deg, q x = 3.7(1), q 2 = 2.1(1), 
r b = 18(1) r g for Case 1; and (ii) E blmd = 6.67(l)keV, O = 
26.7(7) deg, q x = 5.3(1), q 2 = 2.8(1), r b = 4.9(2) r g for Case 3. 
The errors in brackets are evaluated as the 90% confidence re- 
gion for a single interesting parameter when the values of the 
other parameters are fixed. The combination of three parame- 
ters qi, q2, r b thus adjusts the best fit in XSPEC. Nonetheless, 
the clear differences between the models occur consistently with 
theoretical expectations: for Case 3 the lower values of the spin, 
a < 0.87, produce larger^- 2 and the best-fit spin can reach the 
extreme value within the 90% confidence threshold. 
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Of = 30° 
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B f = 60° 


Of = 

f = 30° 


0.9982 
Of = 60° 


1 


56+ 004 


69+ 003 
u - oy -o.04 


QO+0.03 
U - yZ -0.03 


n no 1+0.013 
U ' y51 -0.031 
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66 +0 05 


70 +0 02 

"■'"-0.04 


> 0.966 


> 0.986 


3 


74 +0 05 

' MJ.03 


70 +0 03 

"• ' -0.03 


> 0.991 


> 0.993 



Table 1. The best-fit spin values inferred for the three cases of 
the limb darkening/brightening law, eq. ((H), for the KYRLINE 
model. The artificial data were generated using the KYRLINE 
model with isotropic directionality and the fiducial values of pa- 
rameters (denoted by the subscript "f"). See the main text for 
details. 



3.3. Analysis of simulated data for next generation X-ray 
missions 

In order to evaluate the feasibility of determining the spin of a 
rotating black hole and to assess the expected constraints from 
future X-ray data, we produced a set of artificial spectra. We used 
a simple model prescription and preliminary response matrices 
for the International X-ray Observatory (IXO) missionQ Here 
we limit the energy band in the range 2.5-10 keV. The adopted 
model consists of POWERLAW for continuum (photon index F 
and the corresp onding normalisation Kr), plus KYRLINE model 
dDovciak et al.ll2004bl) for the broad line. The normalisation fac- 
tors of the model were chosen in such a way that the model flux 
matches the flux of MCG-6-30-15. In this section, the simulated 
flux is around 3.1 x 10" 11 ergctrT 2 s _1 with about 3% of the flux 
linked to the broad iron line component. The simulated exposure 
time was 100 ks. 

We generated a set of "fake" spectra (i.e., artificial spectra 
in the XSPEC terminology). These spectra were produced in a 
grid of angular momentum values while assuming isotropic di- 
rectionality, Case 2 in eq. ([5J. We call the assumed angular mo- 
ment of the black hole the fiducial spin and we denote it as «f. 
We performed the fitting loop to these data points using each of 
the three angular emissivity profiles. Once the fit reached conver- 
gence, we recorded the inferred spin a. Figures l9TfT2l show the 
results in terms of best-fit y 2 profiles and the confidence con- 
tours for two different fiducial values of the spin (we assumed 
«f = 0.7, and 0.998). We summarise the values for the inferred 
spin for two inclination angles i = 30deg and i = 60deg in 
Table [TJ 

The fitting procedure was performed in two different ways 
- having the rest of the parameters free or keeping them frozen. 
Obviously the former approach results in an extremely compli- 
cated x 2 space. Therefore, for simplicity of the graphical repre- 
sentation, we plot only the results of the second approach which, 
however, gives broadly consistent results (though it misses some 
local minima of x 1 )- In other words, the plots have the parame- 
ters of the power law continuum, the energy of the line and the 
radial dependence parameter fixed at F = 1 .9, Eq = 6.4 keV, and 
9 = 3. 



4 We used the current version of provisional response matrices 
available at http://ixo.gsfc.nasa.gov/science/responseMatrices.html 
for the glass core calorimeter, dated October 30, 2008. We used the 
energy resolution of 5 eV per bin. 
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Table 2. The best-fit spin and inclination angle values inferred 
for the three cases of the limb darkening/brightening law in the 
KYL3CR model. Data were generated using the KYL2CR model. 
See the main text for details. The quoted errors correspond to the 
90% confidence level. 



The conclusion from this analysis is that the determination 
of a indeed seems to be sensitive within certain limits to the as- 
sumed directionality of the intrinsic emission. The suppression 
of the flux of the reflection component at high values of 6 e may 
lead to overestimating the spin, and vice versa. The middle pan- 
els of Figs. l9TfT2l show the fit results for isotropic directionality, 
which was also the seed model used to generate the test data, 
and so these contours illustrate the magnitude of combined dis- 
persion due to the simulated noise and the degeneracy between 
the spin and the inclination.. The fiducial values are well inside 
the lcr confidence contour in all the graphs in the middle panels. 

However, systematically lower values of the angular momen- 
tum are obtained for the limb brightening profile and, vice versa, 
higher values are found for the limb darkening profile. The mag- 
nitude of the difference is larger for higher values of angular 
momentum. 

3.4. The angular emission profile of the detailed 
reprocessing model 

The results presented in the previous sections show that using 
different emission directionality approaches leads to a differ- 
ent location of the x 1 minimum in the parameter space. Doubt 
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Fig. 9. Test results of the theoretical fits with Of = 0.7, Gf = 30 deg and three different profiles of the emission directionality - 
left: limb brightening (Case 1 in eq. ©), middle: isotropic (Case 2), right: limb darkening (Case 3). The simulated data were 
generated using the POWERLAW + KYRLINE model with isotropic directionality. Top: dependence of the best fit^ 2 values on the 
spin value. The horizontal (dashed) line represents the 90% confidence level. Bottom: contour graphs of a versus O . The contour 
lines correspond to 1,2, and 3 sigma. The position of the minimal value of x 2 is marked with a small cross. The values of x 2 
corresponding to the minimum and to the contour levels are shown at the top of each contour graph. The large cross indicates 
the position of the fiducial values of the angular momentum and the emission angle. Besides the values of spin, inclination and 
normalisation constants of the model, the other parameters were kept fixed at their default values: F = 1.9, £o = 6.4 keV, r m = r ms , 
/"out = 400, q-3, (default values of normalisation constants: Kr = 10~ 2 , Kn ne = 10~ 4 ). 
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Fig. 10. The same as in Figure |9j but for Of = 0.7 and Of = 60 deg. 



about the correct prescription for the emission directionality thus 
brings some non-negligible inaccuracy into the evaluation of the 
model parameters. The magnitude of this error cannot be easily 
assessed as a unique number because other parameters are also 
involved. 



A possible way to tackle the problem is to derive the intrinsic 
spectrum from self-consistent numerical computations. This has 
the potential of removing the uncertainty about the emission di- 
rectionality (although, to a certain degree this uncertainty is only 
moved to a different level of the underlying model assumptions). 
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Fig. 12. The same as in Figure [9J but for Of = 0.998 and Of = 60 deg. 



In this section, we present such results from modelling the artifi- 
cial data generated by numerical simulations, i.e. independently 
of an analytical approximation of the emission directionality pre- 
sented in the previous sections. 

Let us first remind the reader that the orbital speed within the 
inner < 10r g reaches a considerable fraction of the speed of light 
(Fig-ID- Beaming, aberration, and the light-bending all affect the 
emitted photons very significantly in this region. Less energetic 
photons come from the outer parts where the motion slows down 
and the relativistic effects are of diminished importance. This 
reasoning suggests that the analysis of the previous section may 
be inaccurate because the adopted analytical approximations © 
neglect any dependence on energy and distance. 



We applied th e Monte-Carlo radi ative transfer code NOAR 
(see Section 5 in lDumont et alj|2000l) for the case of "cold" re- 
flection, i.e. for neutral or weakly-ionised matter. The NOAR 
code computes absorption cross sections in each layer. Free- 
free absorption and the recombination continua of hydrogen- and 
helium- like ions are taken into account, as well as the direct and 
inverse Compton scattering. The NOAR code enables us to ob- 
tain the angle-dependent intensity for the reprocessed emission. 
The cold reflection case serves as a reference point that we will 
later, in a follow-up paper, compare against the models involving 
stronger irradiation and higher ionisation of the disc medium. 

The directional distribution of the intrinsic emissivity of the 
reprocessing model is shown in Fig.|4](right panel). The contin- 
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Fig. 14. The same plots of data/model ratio, but for a = 0.998 (left: Of = 30°, right: Of = 60°). 



uum photon index F — 1 .9 is considered and the energies are 
integrated over the 2-100 keV range. Although the results of the 
radiation transfer computations do show the limb-brightening ef- 
fect, it is a rather mild one, and not as strong as the Case 1 . In 
the same plot we also show the angular profile of the emissivity 
distribution in different energy ranges: (i) 2-6 keV (i.e. below 
the iron Ka line rest energy); (ii) 6-15 keV (including the iron 
Ka line); (iii) 15-100 keV (including the 'Compton hump'). We 



notice that the energy integrated profile is dominated by the con- 
tribution from the Compton hump, where much of the emerging 
flux originates. However, all of the energy sub-ranges indicate 
the limb-brightening effect, albeit with slightly different promi- 
nence. 

We implemented the numerical results of NOAR modelling 
of a reflected radiation from a cold disc as KYL2CR in the KY 
collection of models. Furthermore, we produced an averaged 
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model, KYL3CR, by integrating KYL2CR tables over all angles. 
Therefore, KYL3CR lacks information about the detailed angu- 
lar distribution of the intrinsic local emission from the disc sur- 
face. On the other hand, it has the advantage of increased com- 
putational speed and the results are adequate if the emission 
is locally isotropic. Furthermore, KYL3CR can be a posteriori 



equipped with an analytical prescription for the angular depen- 
dence (Cases 1-3 in eq.[8]l, which brings the angular resolution 
back into consideration. This approach allows us to switch be- 
tween the three prescriptions for comparison and rapid evalua- 
tion. 
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In order to constrain the feasibility of the aforementioned ap- 
proaches, we generated the artificial data using the POWERLAW 
+ KYL2CR model. The parameters of the model are: photon in- 
dex of the power law T = 1.9 and its normalisation Kr = 0.01, 
spin of the black hole a, inclination angle O , the inner and outer 
radii of the disc r; n = r ms and r out = 400, the index of the ra- 
dial dependence of the emissivity q — 3, and the normalisation 
of the reflection component K\^\2 a =0.1. We simulated the data 
for two different values of the spin, a = 0.7 and a = 0.998, 
and for inclination angles o = 30° and O = 60°. The simu- 
lated flux of the primary power law component is the same as 
in the previous section. However, now an important fraction of 
the primary radiation is reflected from the disc. The total flux 



depends on the extension of the disc and its inclination. Its value 
is 4.9 - 5.6 x 10 _1 1 erg ctrT 2 s _1 for our choice of the parameters. 

As a next step, we replaced KYL2CR by KYL3CR and 
searched back for the best-fit results using the latter model. In 
this way, using KYL3CR we obtained the values of the spin and 
the inclination angle for different directionalities. The fitting re- 
sults are summarised in Table [2] Besides the spin and the in- 
clination angle, only normalisation of the reflection component 
was allowed to vary during the fitting procedure. The remaining 
parameters of the model were kept frozen at their default values. 

The resulting data/model ratios are shown in Figures [T3l and 
[T4l For a = 0.7 and O = 30° the graphs look very similar in 
all three cases. However, the inferred spin value differs from the 
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fiducial value with which the test data were originally created. 
The dependence of the best-fit x 2 statistic on the spin and the 
corresponding graphs of the confidence contours for spin ver- 
sus inclination angle are shown in Figures [TBifLSl again for the 
three cases of angular directionality. These figures confirm that 
for the limb-brightening profile the inferred spin value comes 
out somewhat lower than the correct value, whereas it is higher 
if the limb-darkening profile is assumed. 

In each of the three cases the error of the resulting a de- 
termination depends on the inclination angle and the spin it- 
self. However, we find that the isotropic directionality repro- 
duces our data to the best precision. The limb darkening profile 
is not accurate at higher values of the spin, such as a = 0.998, 
when the resulting x 2 l v value even exceeds 2. The limb dark- 
ening profile is characterised by an enhanced blue peak of the 
line while the height of the red peak is reduced (see Figs. [5]- 
[6]). Consequently, the model profile is too steep to fit the data. 
This is clearly visible in the data/model ratio plots for a = 0.998 
shown in Figures [14] The flux is underestimated by the model 
below a mean energy value E mew of the line (for a = 0.998 and 
i = 30deg E mesLn ~ 5keV) and overestimated above E melLn . This 
fact leads to a noticeable jump in the data/model ratio plot. Only 
the q parameter can mimic this limitation and suit the data by 
increasing its value. It could be the case of the analysis of the 
MCG-6-30-15 data in the previous section. 

A noteworthy result appears in comparing of the contours 
produced by the model with limb brightening and limb dark- 
ening for a = 0.7 and i = 60deg (Fig. [Tol l. Although the for- 
mer model (limb brightening) gives a statistically worse fit with 
X 2 lv = 1.87 than the limb darkening case (x 2 /v = 1.48), the 
inferred values of the spin and the inclination angle are consis- 
tent with the fiducial values within the 3<x level. On the other 
hand, the spin value inferred from the limb darkening model is 
far from the fiducial (i.e., the correct) value. 

4. Discussion 

We investigated whether the spin measurements of accreting 
black holes are affected by the uncertainty of the angular emis- 
sivity law, M(p. e , r e , E e ), in the relativistically broadened iron 
Ka line models. We employed three different approximations of 
the angular profiles, representing limb brightening, isotropic and 
limb darkening emission profiles. For the radius-integrated line 
profile of the disc emission, and especially for higher values of 
the spin, the broadened line has a triangle-like profile. The differ- 
ences among the considered profiles concern mainly the width of 
the line's red wing. However, the height of the individual peaks, 
also affected by the emission directionality, is important for the 
case of an orbiting spot (or a narrow ring), which produce a char- 
acteristic double-horn profile. 

We reanalysed an XMM-Newton observation of MCG-6-30- 
15 to study the emission directionality effect on the broad iron 
line, as measured by current X-ray instruments. We showed the 
graphs of x 2 values as a function of spin for different cases of 
directionality. We can conclude that the limb darkening law fa- 
vors higher values of spin and/or steeper radial dependence of 
the line emissivity; vice versa for the limb brightening profile. 
Both effects are comprehensible after examining the left panel 
of Fig. [5] The limb darkening profile exhibits a deficit of flux in 
the red wing compared with the limb brightening profile. Both 
higher spin value and steeper radial profile of the intensity can 
compensate for this deficit. 

The higher spin value has the effect of shifting the inferred 
position of the marginally stable orbit (ISCO) closer to the black 



hole, in accordance with eq. (@J. Consequently, the accretion disc 
is extended closer to the black hole. The radiation comes from 
shorter radii and it is affected by the extreme gravitational red- 
shift. Hence, the contribution to the red wing of the total disc line 
profile is enhanced. Naturally, these considerations are based on 
the assumption that the inner edge of the line emitting region co- 
incides with the ISCO. Likewise, the steeper radial dependence 
of the emissivity means that more radiation comes from the inner 
parts of the accretion disc than from the outer parts, and this pro- 
duces a similar effect to decreasing the inner edge radius. With 
the limb brightening law the above-given considerations work 
the other way around. 

We further simulated the data with the POWERLAW + KYR- 
LINE model. The simple model allows us to keep better control 
over the parameters and to evaluate the differences in the spin 
determination. We used one of the preliminary response matri- 
ces for the IXO mission and we chose the flux at a level similar 
to bright Seyfert 1 g alaxies observable by current X-ray satellites 
dNandra et al.ll2007t) . The simulations with the KYRLINE model 
confirm that the measurements would overestimate the spin for 
the limb darkening profile and, vice versa, they tend to the lower 
spin values for the limb brightening profile. 

Although the interdependence of the model parameters is es- 
sential and it is not possible to give the result of our analysis in 
terms of a single number, one can very roughly estimate that the 
uncertainties in the angular distribution of the disc emission will 
produce an uncertainty of the inferred inner disc radius of about 
20% for the high quality data by IXO0 We consider this value 
as realistic. 

In the next step, we applied the NOAR radiation transfer 
code to achieve a self-consistent simulation of the outgoing spec- 
trum without imposing an ad hoc formula for the emission an- 
gular distribution. In this paper we assumed a cold isotropi- 
cally illuminated disc with a constant density atmosphere. We 
created new models for the KY suite, KYL2CR and KYL3CR. 
The results of NOAR computations of the cold disc are imple- 
mented in the KYL2CR model, while the KYL3CR model uses 
the angle-integrated tables over the entire range of emission di- 
rections. This enables us to include, a posteriori, the analytical 
formulae for directionality and check how precisely they repro- 
duce the original angle-resolved calculations. We simulated the 
data using the preliminary IXO response matrix and we analysed 
them using the POWERLAW + KYL2CR model. Then we used the 
KYL3CR model to test the analytical directionality approaches 
on these artificial data. 

We found that, for a = 0.998, none of the three assumed 
cases of the directionality profile covered the fiducial values for 
the spin and the inclination angle within the 3cr contour line 
(see Figs. [TTl - fLSl ). The suitability of the particular directional- 
ity prescription depends on the fiducial values of the spin and 
the inclination angle. The limb brightening profile successfully 



5 Here, we refer to the inner disc radius instead of the spin because 
the latter is not as uniform a quantity as the corresponding ISCO radius 
to express the uncertainty simply as a percentage value. However, we 
still suppose that the inner edge of the line emitting region of the disc 
coincides with the marginally stable orbit. In realistic circumstances this 
assumption cannot be satisfied precisely, so it introduc es an additional 
sourc e of error. This has been extensively discussed iBeckwi th et al.1 
2008). The magni tude of the resulting erro r was constrained most re- 
cently in work by Reynolds & Fabian (2008) who applied physical ar- 
guments about the emission properties of the inner flow. It is very likely 
that this discussion will have to continue for some time until the emis- 
sion properties of the general relativistic MHD flows are fully under- 
stood. 
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minimises the x 2 values for a = 0.998 and i = 30 deg, but for 
a = 0.7 and i = 60 deg it gives the worst fit of all studied cases 
of the directionality laws. 

On the whole, we found that the isotropic angular depen- 
dence of the emission intensity fits best. Especially for higher 
values of the black hole spin, the model with the limb darkening 
profile was not able to reproduce the data: the best fit x 2 / v value 
exceeds 5 for a = 0.998 and i = 30 deg, which means a more 
than three times worse fit than using isotropic or limb brighten- 
ing directionality. The inclination angle was underestimated by 
more than 1 deg. 

This is an important result because much of the recent work 
on the iron lines, both in AGN and black hole binaries, has re- 
vealed a significant relativistic br oadening near rapidly rotating 
central black holes ( Milled [2007b . In some of these works, the 
limb darkening law was employed and different options were 
not tested. The modelled broad lines are typically characterised 
by a steep power law in the radial part of the intensity across the 
inner region of the accretion disc, as in the mentioned MCG- 
6-30-15 observation. This behaviour has been interpreted as a 
case of a highly spinning compact source where th e black hole 
rotati onal energy is electromagnetically extracted (Wilm s et al.l 
1200 lh . We conclude that the significant steepness of the radial 
part of the intensity also persists in our analysis, however, the 
exact values depend partly on the assumed angular distribution 
of the emissivity of the reflected radiation. 

It should be noted that, in reality, the angular distribution 
of the disc emission is significantly influenced by the vertical 
structure of the accretion disc. However, our comprehension of 
accretion disc physics is still evolving. In recent years, several 
detailed models have been developed for irradiated black hole 
accretion discs in hydrostat i c equilibrium (see e. g. [Navakshin 
l2000tlBallantvne et alj200"ll lR6zanska et alj|2002l) . These mod- 
els combine radiative transfer simulations with calculations of 
the hydrostatic balance in the stratified disc medium. Aside from 
reprocessing spectra, the models provide solutions for the verti- 
cal disc profile of the density, te mperature, and ionisation frac- 
tions. In Goosmann et al.l (120071) the effects of general relativity 
and adv ection on th e disc m edi um were added. 

In iNavakshinl d2000l) , iGoosmann et all d2007l) . and 
iRozanska & M adei (2008) the reprocessed spectra are evaluated 
at different local emission angles. The shape and normalisation 
of these spectra depend on various model assumptions. Until 
we know in more detail how accretion discs work, it is hard 
to choose which is the "correct" reprocessing model. One 
could always argue that more physical processes and properties 
should be included in the radiative transfer simulations, such 
as the impact of magnetic fields, macroscopic turbulence, or 
different chemical compositions of the medium. Using the 
above-mentioned models for an accretion disc in hydrostatic 
equilibrium may then easily become computationally intense. 

For the practical purpose of data analysis, however, the com- 
putation of large model grids is necessary, which requires suffi- 
ciently fast methods. This is why simple constant density mod- 
els are most often used to analyse the observational data. In fact, 
Balla ntvne et alJ (1200 lh have shown that their reprocessing spec- 
tra for a stratified disc medium in hydrostatic equilibrium can be 
satisfactorily represented by spectra that are computed for irradi- 
ated constant density slabs. Therefore, we include in our present 
analysis the angular emissivity obtained from the modelling of 
neutral reprocessing in a constant density slab. For the purpose 
of this paper we have not discussed in any further detail the de- 
pendence on the ionisation parameter, which we expect to be 
rather important. It will be addressed in a future article. 



We emphasise that the main strategy of the present paper is 
not intended to find the "correct" angular emissivity, as this is 
still beyond our computational abilities and understanding of all 
the physical processes shaping the accretion flow. Instead, we 
examined the three different prototypical dependencies which 
are mutually disparate (i.e., the limb-darkening, isotropic, and 
limb-brightening cases), applied in current data analysis, and 
which presumably reflect the range of possibilities. By includ- 
ing these different cases we mimic various uncertainties, such as 
those in the vertical stratification, and we estimate the expected 
error that these uncertainties can produce in the spin determi- 
nations. Further detailed computations of reprocessing models 
and the angular emissivity are needed in the future in order to 
understand its role in different spectral states of accreting black 
holes. 

5. Conclusions 

Black hole spin measurements using X-ray spectroscopy of rel- 
ativistically broadened lines depend on the definition employed 
of the angular distribution of the disc emission. We studied three 
cases of the directionality profile - limb brightening, isotropic 
and limb darkening. 

We found that isotropic directionality is consistent with our 
radiation transfer computations of the reflection spectra of an 
irradiated cold disc. This eliminates any need for ad hoc ap- 
proximations to the limb-darkening profile. Using an improper 
directionality profile could impact on the other parameters in- 
ferred for the relativistic broad line model. Especially for the 
often used case of limb darkening, the radial steepness can inter- 
fere with the line parameters of the best-fit model by enhancing 
the red wing of the line. 

Uncertainties in the precise position of the inner edge can 
further increase the error of the spin determination; however, it 
appears that the expected magnitude of these errors does not pre- 
vent us from setting interesting and realistic constraints on the 
spin parameter. Our present treatment of the problem is incom- 
plete by neglecting the magnetohydrodynamical effects and their 
influence on the ISCO location. Future improvements in our the- 
oretical understanding of the inner edge location are highly de- 
sirable and will help to improve the confidence in the determi- 
nation of the model parameters. 

We verified our conclusions by comparing them with the re- 
sults based on an XMM-Newton long observation of MCG-6- 
30-15, and we found them to be broadly consistent. The expected 
improvement of sensitivity of the future IXO mission will signif- 
icantly advance the reliability of the spin determinations. 
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